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Abstract 

Imaging the effects of an impurity like Zn in high-T c superconductors [see, e.g., 
S. H. Pan et al., Nature 61 (2000) 746] has rekindled interest in defect problems 
in the superconducting phase. This has prompted us here to re-examine the early 
work of March and Murray [Phys. Rev. 120 (1960) 830] on the linear response 
function in an initially translationally invariant Fermi gas. In particular, we present 
corresponding results for a superconductor at zero temperature, both in the s- and 
in the (i-wave case, and mention their direct physical relevance in the case when the 
impurity potential is highly localized. 
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1 Introduction 

The unconventional nature of the pairing state in the high-T c superconductors 
(HTCS) has been recently provided with further evidence by the observation of 
the spatial inhomogeneity of the electron distribution around an isolated non- 
magnetic impurity, such as Zn, in cuprate Bi-2212, imaged by scanning tun- 
neling microscopy (STM) [1,2,3,4]. These experimental studies demonstrated 
the existence of impurity-induced quasi-bound states (resonances) near the 
Fermi level. Such results have been naturally interpreted as fingerprints of the 
momentum-space anisotropy of the order parameter in the HTSC In partic- 
ular, the four-lobed structure of the STM differential conductivity observed 
around an isolated impurity is a clear manifestation of the (i-wave symmetry 
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of the superconducting gap in these materials. More recently, similar results 
have been reported for Nd/Ba-substituted YBCO thin films [5,6], thus lending 
further support to the ci-wave scenario for the hole-doped high-T c cuprates. 

The idea that an anisotropic superconducting gap should give rise to directly 
observable spatial features in the tunneling conductance near an impurity was 
suggested by Byers et al. [7], whereas earlier studies [8] had considered per- 
turbations of the order parameter in unconventional superconductors to occur 
within a distance of the order of the coherence length £ around an impurity. 
Later, it was shown that an isolated impurity in a rf-wave superconductor pro- 
duces virtual bound states close to the Fermi level, in the nearly unitary limit 
[9]. Such a quasi-bound state should appear as a pronounced peak near the 
Fermi level in the local density of states (LDOS) corresponding to an impurity 
site [10], as is indeed observed in Bi-2212 [1] and YBCO [5]. 

The location and number of such low-energy peaks in the LDOS around an 
impurity could also provide important new insights into the nature of the 
pseudogap state in the underdoped cuprates above the critical temperature T c 
[11,12]. Within the preformed pairs scenario, a pseudogap in the normal state 
is related to the presence of fluctuating pairs, still lacking phase coherence, 
characterized by a much shorter coherence length and lifetime than supercon- 
ducting pairs below T c [13]. On the other hand, the ci-density wave (DDW) 
scenario has been recently proposed [14], in which the normal state is charac- 
terized by an actual gap in the quasiparticle spectrum, due to a hidden broken 
symmetry of d x 2_ y 2-tjpe. In the presence of a pseudogap depleting the DOS 
at the Fermi level, an impurity-induced resonant state should survive above 
T c , although broadened in energy. However, if preformed pairs above T c have 
a superconducting origin, they would be composed of a superposition of hole 
as well as electron states, like Bogoliubov quasiparticles below T c . Therefore, 
impurity-induced resonances in the LDOS probed by STM spectroscopy would 
be expected at both signs of the applied bias. On the other hand, if the normal 
state (pseudo)gap has no direct connection with the superconducting g ap, as 
in the DDW scenario, then only a definite kind of peaks, either of electron or 
of hole kind, should appear. 

The spatial inhomogeneity of the local electronic structure around an isolated 
impurity in an unconventional superconductor has been already studied via 
the solution of the Bogoliubov-De Gennes equations for a disordered system 
[15,16,17,18], and within a self-consistent approach, taking into account the 
changes in the order parameter induced by the impurity potential [19]. Here, 
we follow a different approach, and address the problem of finding the change 
in the local electron density around an isolated impurity in an anisotropic 
superconductor, within linear response theory. 

The outline of the paper is then as follows. After reviewing the expression of 
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the linear response function F in real space for a uniform electron gas [20], in 
Sect. 2 we calculate the first-order Dirac density matrix, both for an isotropic, 
s-wave superconductor in three dimensions (3D), and for an anisotropic, d- 
wave superconductor in two dimensions (2D). In Sect. 3 we then generalize the 
above expression for F to the superconducting case, but now in momentum 
space, where it is directly related to the Fourier transform of the electron 
density change in a superconductor around a highly localized impurity. A 
numerical analysis then demonstrates that a rf-wave order parameter gives 
rise to an azymuthal modulation in the momentum dependence of F, which is 
responsible of the four-lobed pattern observed in the density change around an 
impurity. Later in Sect. 4 we summarize and give directions for future work. 



2 Linear response theory for normal and superconducting state 



The linear response function F(r, r', Ep) for a one-body potential V(r) intro- 
duced into an initially uniform Fermi gas provides the density change Sp(r) 
in that gas due to the 'impurity' generating V(r). This quantity F is given in 
the early work of March and Murray [20] and is translationally invariant and 
spherically symmetric, namely 



where ji(x) = (sin x — x cos x) jx 1 is the spherical Bessel function of order one, 
kp the modulus of the Fermi momentum, and Ep = \kp the Fermi energy (in 
units such that h — m — 1, to be used throughout the present work). Here, 
prompted by the rekindling of interest in defect problems in superconductors 
[1,2,3,4], the analogous response function is calculated within the BCS model 
of the superconducting phase [21], for both isotropic and anisotropic super- 
conductors. Due to the 'rounding off' of the Fermi momentum distribution 
n(p), which at T = in the normal state is the usual step function terminat- 
ing at Fermi momentum pp = hkp, the superconducting (S) response function 
Fs(r — r', Ep) has more pronounced 'damped' oscillations at large separations 
|r — r'|, as compared to the normal state case, the damping factor of the os- 
cillations measuring the 'blurring' of the 'edge' at n{pp) due to Cooper pair 
formation. In addition to that, for an anisotropic superconducting state, as is 
the case for a d-wave order parameter, the superconducting response function 
F s (r — r', Ep) is expected to lose spherical symmetry. 

Returning to the study of March and Murray [20] , these authors obtained the 
displaced charge Sp(r) due to a perturbation V(r) introduced into an initially 
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uniform Fermi gas having electron density po = k F /?>Tx 2 as 

Sp(r) = J V(r')F(\r-r'\,E F )d 3 r'. 



(2) 



Here, our interest is to generalize Eq. (2) in order to treat the effect of a given 
perturbation V(r) introduced into the superconducting phase as described by 
the BCS model [21]. Then we can write the analogue of Eq. (2) in the form 



5p s (r) = J V(r')F s (v-r',E F )dh'. 



(3) 



To proceed to calculate F s below, we note next the result of Stoddart, March, 
and Stott [22] (see also Ref. [23]) that 



dF(r,r',E) 
dE 



= 2 Re 



G(t,t',E + ) 



fry(r,r',£) 
dE 



(4) 



where E + = E + iS. Inserting the Green function for an outgoing spherical 

wave 



G(r, r', E + ) 



exp(i/c|r — r'|) 
47r|r — r'l 



(5) 



where E = \k 2 , and the first-order Dirac density matrix 



7 (r, r', E F ) = 2 £ exp[ik • (r - r')] = J 1 (6) 

|k|<fc F 



k F ji(k F 


|r-r'|) 


7T 2 k F \ 


r — r' 





into Eq. (4), one regains the response function, Eq. (1), after some manipu- 
lation. We also remind that Eq. (4) can be equivalently cast in the integral 
form [24]: 

F(r, r', E F )=iJ G(r, r', E)G(r', r, E) dE. (7) 



In the superconducting state, the upper branch of the Bogoliubov excitations 
is characterized by a dispersion relation E k = + A£, where £ k = \{k 2 — k F ) 
is the free particle dispersion, measured with respect to the Fermi level, and 
Ak is the superconducting gap at T = (in general, a function of momentum 
k). Eq. (6) is immediately generalized then to read in the superconducting 
phase 

75 (r, r', E F ) = 2 ]T n s (k) exp[ik • (r - r')], (8) 

k 
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Fig. 1. Shows ratio of first-order density matrix to its diagonal value in 3D for three 
different values of the ratio Aq/E f = 0.1, 0.3, 0.5. Solid line is normal state result. It 
should be actually emphasized that Aq/E f < 1 for most cuprates. The larger values 
of the parameter Ao / E F used here and in Fig. 2 below are intended to magnify the 
effect of a superconducting gap in the first-order Dirac density matrix. 



where ng(k) denotes the momentum distribution in the superconducting phase 
at T = 0, which coincides with the coherence factor t£ = |(1 — E k ) of 
BCS theory [21]. In the s-wave case, corresponding to having a constant gap 
Ak = A over momentum space, is characterized by an inflection point 
at k F of width ~ 2A (see Fig. 2.4 in Ref. [21], or Fig. 7.12 in Ref. [25]), 
and 7s(r, r',E F ) = 7s(|r — r'\,E F ). We have plotted Eq. (8) for 7s(|r — 
r'|, £'f)/7s(0, E f ) as a function of separation |r — r'| for three instances of 
the ratio A /E F in Fig. 1, where these three curves are compared with the 
normal state result, Eq. (6). The 'damping' of the oscillations in the normal 
state response function on going to the superconducting phase is evident from 
the plots in Fig. 1. 

The above derivation can be readily generalized to an anisotropic super- 
conductor. In the case of the high-T c cuprates, we assume a d-wave gap 
Ak = A cos2#, with k now a two-dimensional (2D) wavevector, 9 being 
the angle it forms with the x direction. Such an order parameter is character- 
ized by nodes for 6 = ±7r/2, corresponding to the directions y = ±x. In the 
normal case, but now in 2D, we find the translationally invariant, spherically 
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symmetric first-order Dirac density matrix: 

,(ry,E F) = k l J ff-^ , (9) 
7r k F \r — r\ 

with Ji(x) here denoting the Bessel function of the first kind and order one, 
which again displays oscillations as Eq. (6), but now with a weaker asymp- 
totic decay in the limit kp\r — r'| ^> 1 (see Fig. 2, upper left panel). In the 
superconducting case, taking into account the ci-wave symmetry of the or- 
der parameter, 7s(r — r',Ep), Eq. (8), loses its spherical symmetry. We plot 
7s(r — r', E F )/^s(0, E F ) in Fig. 2 as a function of the 2D vector r — r'. As in 
the isotropic, 3D case, we recover a 'damping' of the oscillations in the normal 
state response function in going into the superconducting phase, which is now 
more pronounced in the fully gapped x and y directions, than in the gapless 
directions y = ±x. 



3 Linear response function for an anisotropic superconductor 



We now turn to the problem of explicitly evaluating the linear response func- 
tion for a 2D superconductor, both in the s- and in the ci-wave case. We start by 
noting that Eqs. (2-3) express the relation between the density change bp and 
the impurity potential V as a convolution in real space, the convolution kernel 
being the linear response function F. Such a convolution naturally translates 
into a simple product in momentum space, 5p(q) = V(q}F s (q, E F ). Moreover, 
in the case of a highly localized impurity, as is the case of a 5-function-like 
impurity potential in real space [7], one has V(q) = Vq, so that the Fourier 
transform of the linear response function -Fs(q, Ep) practically coincides with 
the Fourier transform of the density change. In the normal state, the latter 
quantity is readily obtained by Fourier transforming Eq. (7) as a convolution 
of Green's functions in momentum space. In the superconducting state, such 
an expression is generalized by the static limit (u = 0) of the polarization 
function, which for a 2D superconductor reads [26] 



i2l 

F 5 (q,£ F ) = ^-Tr JdEj -— pT 3 G(k,£)T 3 G(k- q ,£) 

= 9 f K+q^k + ^k+q) 2 nm 

J (27T) 2 ^ £ k+ q + £k ' ^ ) 

where again and u\ = 1 — are the coherence factors of BCS theory 
[21], but now in general depending on an anisotropic gap Ak- In Eq. (10), G 
denotes the matrix Green's function in Nambu notation [21], and r 3 is a Pauli 
matrix. 
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Fig. 2. Shows ratio of first-order density matrix to its diagonal value for three 
different values of the ratio Aq/Ey = (normal case), 0.1, 0.3, 0.5 (left to right, 
bottom to top), but now for the 2D, d-wave case. As in Fig. 1, all contour plots are 
against k F (r — r') (now a 2D vector). Damping of the normal-state oscillations is 
more pronounced in the fully gapped x and y directions than in the gapless y = ±x 
directions. 

In the normal state (A — > 0), the integration in Eq. (10) is feasible analyti- 
cally, and we find: 

( -, if q < 2k F , 

F(q, E F ) = F 5 (q, E F ; A = 0) = f 7-^4 (11) 

, if q > 2k F . 

\ TT q 

Therefore, we find that the linear response function for a 2D system in the 
normal state is characterized by a step discontinuity at q — 2k F (Fig. 3), which 
is responsible of damped Friedel-like oscillations of period ~ Tr/k F for F(r, E F ) 
in real space. 

In the superconducting state, for an s-wave, isotropic order parameter in 2D, 
such a discontinuity is smoothed over a width of ~ 2A around q = 2k F , as 
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Fig. 3. Linear response function in momentum space, F$(q, Ep), Eq. (10), in the case 
of a translationally invariant, isotropic 2D system. Dashed line shows the normal 
state analytical result, Eq. (11). Continuous lines show our numerical results for the 
superconducting state, when an s-wave order parameter is assumed, for the same 
values of Aq/Ef as in Fig. 1. Notice the step discontinuity at q = 2kp in the normal 
state, which gets smoothed over a width ~ 2Aq in the superconducting case. 



is shown numerically in Fig. 3. This is again expected to give rise to more 
pronounced damped oscillations in the r-dependence of the linear response 
function in real space. 



In the case of a 2D, rf-wave superconductor, F 5 (q, E F ) must be analyzed as 
a function of q as a vector. Numerical results are plotted in Fig. 4 for dif- 
ferent values of the ratio A /Ep, and compared to the (isotropic) result for 
the normal state. As in the isotropic, s-wave case, the opening of a supercon- 
ducting gap at the Fermi level tends to smear out the step discontinuity at 
q = 2kp over a width ~ 2A . However, such an effect is more enhanced in the 
q x and q y directions, corresponding to maxima in the energy gap, than in the 
q y = ±q x directions (gap nodes), where it is virtually absent. This gives rise 
to an azymuthal modulation of -Fs(q, Ep), of period tt/2, which is ultimately 
responsible for the clover pattern in the density change around an isolated 
impurity, as observed in STM imaging experiments. 
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Fig. 4. Linear response function Fs(q, Ep), Eq. (10), for a d-wave superconductor 
in 2D. Actually displayed are contour plots of Fg(q,Ep) as a function of q/kp for 
Aq/Ef = 0.0 (normal state), 0.1, 0.3, 0.5 (left to right, top to bottom). In going into 
the superconducting state, the step discontinuity at q = 2kp gets smoothed more 
in the q x and q y directions (corresponding to gap maxima) than in the q y = ±q x 
directions (corresponding to gap nodes). 

4 Conclusions and directions for future work 

Motivated by the recent imaging results around an isolated impurity in HTSC 
[1,2,3,4,5,6], we have analyzed the linear response function in momentum space 
for a 2D uniform electron gas, both in the normal and in the superconducting 
state. The opening of an energy gap at the Fermi level manifests itself already 
in the radial dependence of the first-order Dirac density matrix. In the s-wave 
case, its diagonal element 7s(|r — r'\,E F ) decays more rapidly with distance 
from the diagonal than does its normal state counterpart. In the d-wave case, 
fingerprints of the anisotropic momentum dependence of the energy gap are 
already present in the angular dependence of 7s(r — r',Ep), with more pro- 
nounced damping of the oscillations in the fully gapped x and y directions 
than in the gapless y = ±x directions. Such a behaviour is confirmed by 
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the momentum dependence of the linear response function Fg(q, Ep)- In the 
normal, isotropic state, the latter function is characterized by a step discon- 
tinuity at q = 2k F , giving rise to (damped) Friedel-like oscillations in its real 
space counterpart. The discontinuity at q = 2/c F gets smoothed over a width 
~ 2A in going into the superconducting state. For a rf-wave superconduc- 
tor, however, such an effect is more pronounced in the fully gapped q x and q y 
directions, than along the q y = ±q x nodes. This gives rise to an azymuthal 
modulation of Fg(q, E F ), which is ultimately responsible of the four-lobed 
pattern in the density change, observed in STM imaging experiments around 
an impurity. Before a direct comparison could be made with such experimen- 
tal results, however, we believe that a more accurate band dispersion relation 
should be considered (see, e.g., Ref. [27]), as well as a more realistic choice for 
the impurity potential, which we reserve to future work. 



Acknowledgements 

The authors are indebted with M. Salluzzo, F. Siringo, and J. Verbeeck for 
helpful discussions. One of us (N. H. M.) made his contribution to this work 
during a visit to the Physics Department, University of Catania, in the year 
2000. Thanks are due to the Department for the stimulating environment 
and for much hospitality. N. H. M. also wishes to thank Professor V. E. Van 
Doren for his continuing interest and support. G. G. N. A. acknowledges the 
University of Antwerp (RUCA) for much hospitality during the period in which 
this work was revised and brought to completion. 



References 

[1] S. H. Pan, E. W. Hudson, K. M. Lang, H. Eisaki, S. Uchida, J. C. Davis, 
Imaging the effects of individual zinc impurity atoms on superconductivity in 
Bi 2 Sr 2 CaCu 2 8+5 , Nature 403 (2000) 746. 

[2] E. W. Hudson, S. H. Pan, A. K. Gupta, K. W. Ng, J. C. Davis, Atomic-scale 
quasi-particle scattering resonances in Bi 2 Sr 2 CaCu 2 08+<5, Science 285 (1999) 
88. 

[3] E. W. Hudson, K. M. Lang, V. Madhavan, S. H. P. H. Eisaki, S. Uchida, J. C. 
Davis, Interplay of magnetism and high-T c superconductivity at individual Ni 
impurity atoms in Bi 2 Sr 2 CaCu 2 8+5 , Nature 411 (2001) 920. 

[4] S. H. Pan, J. P. O'Neal, R. L. Badzey, C. Chamon, H. Ding, J. R. Engelbrecht, 
Z. Wang, H. Eisaki, S. Uchida, A. K. Guptak, K. W. Ng, E. W. Hudson, 
K. M. Lang, J. C. Davis, Microscopic electronic inhomogeneity in the high- 
Tc superconductor Bi 2 Sr 2 CaCu 2 8+a; , Nature 413 (2001) 282. 



422 



[5] M. Iavarone, M. Salluzzo, R. Di Capua, M. G. Maglione, R. Vaglio, 
G. Karapetrov, W. K. Kwok, G. Crabtree, STM tunneling spectroscopic studies 
of YNd a; Ba2_xCu307_ 5 thin films, submitted to Phys. Rev. B ... (2002) ... 

[6] M. Salluzzo, A. Andreone, F. Palomba, G. Pica, R. Vaglio, I. Maggio- 
Aprile, 0. Fischer, A systematic study of the role of Nd/Ba disorder on the 
superconducting properties of Ndi +:r Ba2-xCu307_a, Eur. Phys. J. B 24 (2001) 
177. 

[7] J. M. Byers, M. E. Flatte, D. J. Scalapino, Influence of gap extrema on the 
tunneling conductance near an impurity in an anisotropic superconductor, Phys. 
Rev. Lett. 71 (1993) 3363. 

[8] E. S. Choi, P. Muzikar, Impurity magnetic moments in unconventional 
superconductors, Phys. Rev. B 41 (1990) 1812. 

[9] A. V. Balatsky, Impurity induced virtual bound states in <i-wave 
superconductors, Phys. Rev. B 51 (1995) 15547. 

[10] M. I. Salkola, A. V. Balatsky, D. J. Scalapino, Theory of scanning tunneling 
microscopy probe of impurity states in a ci-wave superconductor, Phys. Rev. 
Lett. 77 (1996) 1841. 

[11] H. V. Kruis, I. Martin, A. V. Balatsky, Impurity-induced resonant state in a 
pseudogap state of a high-T c superconductor, Phys. Rev. B 64 (2001) 054501. 

[12] Qiang-Hua Wang, Mechanism of pseudogap probed by a local impurity, Phys. 
Rev. Lett. 88 (2002) 057002. 

[13] M. Randeria, Precursor pairing correlations and pseudogaps, in: G. Iadonisi, 
R. J. Schrieffer, M. L. Chiofalo (Eds.), Proceedings of the CXXXVI 
International School of Physics "E. Fermi" on "Models and phenomenology 
for conventional and high-temperature superconductivity", IOS, Amsterdam, 
1999, preprint cond-mat/9710223. 

[14] S. Chakravarty, R. B. Laughlin, D. K. Morr, C. Nayak, Hidden order in the 
cuprates, Phys. Rev. B 63 (2001) 094503. 

[15] A. Ghosal, M. Randeria, N. Trivedi, Spatial inhomogeneities in disordered d- 
wave superconductors, Phys. Rev. B 63 (2000) R020505. 

[16] S. Haas, K. Maki, Quasiparticle bound states around impurities in d x 2_ y 2-wme 
superconductors, Phys. Rev. Lett. 85 (2000) 2172. 

[17] K. Maki, S. Haas, Impurity bound states and symmetry of the superconducting 
order parameter in Sr 2 Ru0 4 , Phys. Rev. B 62 (2000) R11969. 

[18] S. Haas, K. Maki, Temperature dependence of impurity bound states in d x 2_ y 2- 
wave superconductors, submitted to Europhys. Lett. ... (2000) preprint 
cond-mat/0011004. 

[19] M. E. Flatte, J. M. Byers, Local electronic structure of defects in 
superconductors, Phys. Rev. B 56 (1997) 11213. 



423 



[20] N. H. March, A. M. Murray, Relation between Dirac and canonical density 
matrices, with applications to imperfections in metals, Phys. Rev. 120 (1960) 
830. 

[21] J. R. Schrieffer, Theory of Superconductivity, W. A. Benjamin, New York, 1964. 

[22] J. C. Stoddart, N. H. March, M. J. Stott, Electronic structure of point defects 
and impurities in metals, Phys. Rev. 186 (1969) 683. 

[23] W. Jones, N. H. March, Theoretical Solid-State Physics. Non-equilibrium and 
Disorder, Vol. 2, Dover, New York, 1986. 

[24] F. Flores, N. H. March, Y. Ohmura, A. M. Stoneham, Asymptotic displaced 
charge round impurities in metal crystals with and without surfaces, J. Phys. 
Chem. Solids 40 (1979) 531. 

[25] N. H. March, W. H. Young, S. Sampanthar, The Many-Body Problem in 
Quantum Mechanics, Dover, New York, 1995. 

[26] R. E. Prange, Dielectric constant of a superconductor, Phys. Rev. 129 (1963) 
2495. 

[27] I. Martin, A. V. Balatsky, Local structure of impurity states in (i-wave 
superconductors, Physica C 341-348 (2000) 125. 



424 



